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Carrier injection into carbon nanotubes and graphene nanoribbons, contacted by a metal coating 
over an arbitrary length, is studied by various means: Minimal models allow for exact analytic 
solutions which can be transferred to the original system with high precision. Microscopic ah initio 
calculations of the electronic structure at the carbon-metal interface allow us to extract — for Ti and 
Pd as contacting materials — realistic parameters, which are then used in large scale tight-binding 
models for transport calculations. The results are shown to be robust against nonepitaxially grown 
electrodes and general disorder at the interface, as well as various refinements of the model. 
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I. INTRODUCTION 

The high electrical conductivity of metallic car- 
bon nanotubes (CNTs) can be attained, thanks to a 
unique combination of several features. The quasi-one- 
dimensional crystal structure, together with a low den- 
sity of defects, allows us to explore the theoretical limit 
of conductance of Ae'^/h at the charge neutrality point. 
The stiffness of carbon-carbon bonds reduces the effect of 
electron-phonon coupling at room temperature .^ Also, 
restricting electron movement to a single dimension re- 
sults in a very small phase space, which strongly reduces 
the effectiveness of scattering. A further reduction of 
backscattering is caused by the low density of states at 
the Fermi energy in combination with a high Fermi ve- 
locity. Considering all these factors, measured ballistic 
lengths of several microns^ become understandable. Yet, 
to exploit the potential for carrying current densities of 
up to 10^ A/cm^,'' the contacts at a nanometer scale be- 
come crucial. 

For graphene, which has been shown to allow ballis- 
tic transport over similar lengths,''^ the contacts play an 
equally important role.^ Finite-width graphene nanorib- 
bons (GNRs) are discussed as nanoelectronic devices and 
interconnects, showing both similarities and distinctive 
differences when compared to nanotubes. With first 
experimental results becoming available,^ many of the 

theoretical predicions about GNRs will soon be put to 

test7,8,9,io 

Experimentally, a crucial factor for obtaining good 
metallic contacts are the wetting properties of the ma- 
terial. Thus, it has been observed that Ti, Ni and 
Pd form continuous coatings on single-wall CNTs while 
Au, Al and Fe form isolated particles. "'^^ Furthermore, 



among several common contact metals, Ti was found to 
be the only one where true chemical bonds could be ob- 
served, while the others showed only weak van der Waals 
interactions.^^ Surprisingly enough, Pd — traditionally 
known as a rather poor conductor — was found to form 
better and more reliable Ohmic contacts for CNTs than 
Ti^ and could be successfully applied to produce a CNT 
field effect transistor with Ohmic contacts. ^"^ It is gener- 
ally believed that this superiority of Pd is due to its high 
work function ((/)pd = 5.1 eV) that matches well with 
that of CNTs [e.g., 0(7, o) cnt = 5.1 eV (Ref. 14)] and 
thereby avoids a high Schottky barrier. Pt, which has an 
even higher work function, would therefore be expected 
to perform even better as a contacting material, but as 
it turns out it does not form Ohmic contacts at all.^^ 

A further experimental puzzle is the question of the 
effective length of contacts formed by a coating metal 
layer: While some studies report that transport occurs 
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FIG. 1: (Color online) Minimal model for extended contacts 
solvable analytically: a linear chain of identical atoms (one 
orbital per atom) with hopping integral 7 between nearest 
neighbors. To the left, the chain continues infinitely; at the 
right end, TV atoms are contacted, each by an independent 
wideband lead of strength A. For defining the transmission, 
the system is virtually split into three regions: the conductor 
C and the leads L and R. 



only at the edge of the contact,^ others state that the 
contact resistance depends on the length of the contact^*' 
or that extended contacts show a clearly distinct behavior 
from pointlike contacts. ^^ 

Various theoretical studies have been conducted to in- 
vestigate these issues: In ab initio studies comparing 
Au, Pd and Pt contacting an (8, 0) CNT, Pd was found 
to have the lowest Schottky barrier. ^^ Another ab ini- 
tio study comparing the metal-graphene bonding of the 
same three metals indicated a very small binding energy 
for Au. For Pd it is somewhat larger, while for Pt the 
bonding is yet stronger. ^^ The bad contacts formed by Pt 
are here attributed to a clustering effect of larger metal 
grains. Direct ab initio simulations of transport in a 
metal-contacted CNT compared Pd and Au,^°'^^ find- 
ing again that Pd forms superior contacts. Ti was also 
found to form strong bonds to a nanotube surface^^ and 
to form superior point contacts for CNTs.^'^ 

Apart from these practical issues in explaining and im- 
proving the quality of materials, the study of contact 
models is also of great theoretical relevance: In studying 
the physics of electronic devices at the nanometer scale, 
it is generally crucial to have detailed control over the 
contacts. ^'*'^^ Indeed, specifying properties of nanoelcc- 
tronic devices is generally completely meaningless with- 
out clearly stating the way the system was contacted or — 
for theoretical studies — how the contact was modeled. 

In a previous work,^^ we have demonstrated the coun- 
terintuitive finding that the optimal transparency in ex- 
tended contacts for CNTs is achieved not by strong 
chemical bonding — as it would be the case for pointlike 
contacts^"^'^^''^®''^^ — but rather by contact materials that 
gently couple to the tube surface, allowing us to exploit 
the length of the tube-metal interface to smoothly inject 
the electrons with minimal reflection at the contact. 

In this paper, we will present a detailed analysis of the 
model used for our previous findings. A minimal model, 
reducing the CNT or GNR to a single atomic chain and 
considering a semi- infinite wire with only a single contact 
at one end, allows an analytic solution that gives detailed 
insight in the mechanism of an extended contact (see 
Fig. 1). The results from this model have a close relation 
to the physics of Breit-Wigner resonances, which will be 
briefiy sketched out at the beginning. After a thorough 
analysis of the minimal model, we will demonstrate how 
the results can be transferred to more realistic structures 
and explain in detail how the missing parameters could 
be quantiatively estimated from ab initio calculations of 
Pd and Ti as two typical contacting materials. 



rates Fl and Fr by the Breit-Wigner equation' 
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The bell-shaped peak in this expression as a function of 
the energy is well known. What is rarely noted in litera- 
ture, however, is the fact that also for fixed energy E and 
one fixed contact Fl, the conductance as a function of the 
other contact G (Fr) has a bell shape with an optimum 
at the balanced coupling Fr = Fl. A new perspective 
to this old issue was provided in the recent experiments 
by Griiter et al.'^'^ Small couplings (Fr <^ Fl) result in a 
linear Fr dependence of the conductance typical of tun- 
neling phenomena, while for large coupling (Fr ^ Fl), 
such better contact Fr results in an overall suppressed 
conductance. One way to understand this counterintu- 
itive behavior is to consider the tunneling rate F as a 
measure for the chemical bond between the conducting 
orbitals of the molecule and the lead: A strong bond to 
one of the leads causes the molecule itself to virtually be- 
come part of that lead so that we observe the physics of 
a single point contact. Furthermore, the strong bonding 
results in a strong redistribution of the spectral weight 
of the energy level in the molecule, i.e., in a low local 
density of states (LDOS) at the energy e. The tunnel- 
ing conductance, which directly probes this LDOS, will 
therefore be suppressed by large Fr. 

Once the length A'^ of a contact is increased for an ex- 
tended molecule, the optimal value decreases monotoni- 
cally with the number of contact points N, as displayed 
in Fig. 2. As will be shown later [see Eq. (7)] this value 
scales like Fr = Fl In N/N for large iV. 



III. ANALYTICAL MODEL 

A minimal model that captures the essential physics 
of extended contacts is set up as follows. The CNT or 
GNR is represented by a linear chain of atoms with the 
hopping integral 7 and the onsite energy e = (fixing the 
energy offset). A two-probe setup is defined by selecting 
an arbitrary single electron as the "conductor" and the 
semi-infinite sections at both ends as "leads". In this 
unmodified setup, the system is fully transparent, so the 
transmission T (E) is equal to the number of channels 
A'^ch at any given energy. The single cosine-shaped band 
of the linear chain provides a single transmission channel 

n,nd{E) = e (s + 27) e (-£; - 27) 



II. BREIT-WIGNER RESONANCE 

The conductance for a molecular junction, consisting 
of a single energy level e in a two-terminal setup between 
two leads, is given in terms of the left and right tunneling 



which presents a theoretical upper transmission limit 
when scattering at the contacts could be neglected. 

An "extended contact" to the linear chain is now mod- 
eled by replacing the semi-infinite lead by a finite iV-atom 
chain contacted in each atom individually by a wideband 
lead of strength A. (For a sketch of the model, see Fig. 1.) 
A full solution of this model is obtained by calculating 
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FIG. 2: (Color online) Breit-Wigner resonance in the conduc- 
tance of an extended molecule sandwiched between two leads. 
The internal hopping is fixed as 27 = Fl to allow an optimal 
match to the left contact. In the case A'' = 1, the system is 
identical to the molecular junction described by Eq. (1). The 
transmission shows the shift of the Breit-Wigner peak towards 
lower Fr with growing A'^. The functional form of this shift 
can be approximated for large N as Fr = Fr In N/N [see 
Eq. (7)]. 



the conductance as the quantum mechanical transmission 
probabiHty, 

G{E) = '^T{E). 

This can be done within the Landauer approach to trans- 
port by means of the Green function formahsm, as shown 
in the Appendix. As a result, one obtains (for e = 0) 



T{E) = 



S^A-Ey^nm [fN iE/2-f - iA/47)] 



S/7 - 1^4-^2/72 - 2fN (S/27 - iA/47) 



with f]\[ {x) — Un-1 {x) /Un (x). Un{x) are the Cheby- 
shev polynomials of the second kind, as given in Eq. (A4). 
To gain a full understanding of the physics described by 
this expression, the transmission T{E) is plotted as a 
function of the energy for different values of the two pa- 
rameters N and A (i.e., the length and the quality of 
the contact region) in Fig. 3. Two regimes can be iden- 
tified: An N -resonant regime for low A/small N , where 
the transmission shows about as many peaks as there 
are atoms in the contact region, and an N -independent 
regime for high A /large N , where the transmission shows 
no resonances and depends only on A. 

The two different plots in Fig. 3 illustrate two as- 
pects: For fixed N with increasing A, the transparency 
of the system improves, goes through an optimum 
point, and degrades again, while for fixed A, the trans- 
parency improves with growing N and saturates at an 
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FIG. 3: (Color online) Transmission through the system dis- 
played in Fig. 1, as given in Eq. (A8). Top panel: for a fixed 
contact length iV, starting from low A, the transmission first 
improves, reaches an optimum, and then degrades again at 
high A. Bottom panel: For fixed contact strength A, trans- 
mission improves with growing A'' and saturates for large A'". 



iV-independent optimum. In both cases, the transmis- 
sion goes through the two different regimes. 

Especially the last point can be seen more clearly by 
looking at the reflection i? = 1 — T in the energy range 
of the single channel of our system. 



R = 



E + 1^472 -E^- 27/Ar (g/27 - iA/47) 
E - iv/472 - £;2 - 27/Ar (£;/27 - iA/47) 



(2) 



In Fig. 4, this observable is plotted in a logarithmic scale, 
illustrating that the average value of the transmission 
already saturates at A^ = 100 (specifically for A = O.I7). 
For larger values of A^, the overall transparency is not 
improved any further, but the A'^-dependent resonances 
are smoothed out. 




FIG. 4: (Color online) Contact reflection in the system dis- 
played in Fig. 1, as given in Eq. (2). At fixed contact strength 
A, with growing A*', the contact becomes more transparent 
and saturates at an A^-independent value. 



To better understand the origin of this saturation, we 
fix the energy to the half filling case E — e and study 
the transmission and the reflection for varying contact 
lengths N (see Fig. 5): 



R{E^e) 



1 + i/AT (-iA/47) 



1 - iJN (-iA/47) 



(3) 



Similar data were numerically obtained before for 
CNTs.'^'* A different insight, however, can be gained from 
the reflection in a logarithmic scale: Ignoring the even- 
odd oscillations in N , one observes flrst an exponential 
decay of R with increasing iV, followed by an abrupt 
crossover to an A^-independent value. Both the rate of 
decay and the saturation value depend on A in such a 
way that for lower values of A, the transparency initially 
improves more slowly with the contact length, but ulti- 
mately R saturates at a lower value, which means higher 
contact transparency. This result, which we presented 
before based on numerical calculations on CNTs,^^ will 
be studied in more detail in the following using our ana- 
lytical expressions. 

An expression for the iV-independent regime can easily 
be obtained as the limit N ~* 00 oi Eq. (A6) as 
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which can be further simplifled for A <C 7 to obtain 

A2 



(4) 




FIG. 5: (Color online) Transmission T (top panel) and con- 
tact reflection R (bottom panel) in the system displayed in 
Fig. 1 at fixed energy E = e for varying contact length A^ and 
selected values of the contact strength A. 



The validity of this approximation is illustrated in Fig. 6. 
The approach for flnding the corresponding approxi- 
mation for the iV-resonant regime is less rigorous since 
a simple limit is not sufficient to capture the behavior 
in this case. A far better approximation is found graph- 
ically: The straight section in a semilogarithmic scale 
plot (inset of Fig. 6) indicates a clean exponential law. 
The missing coefficients are easily found from a Taylor 
expansion in A = 0, yielding 



resonant 
:N odd 
resonant 



exp (-iVA/7) , 
exp[-(Ar4-l)A/7] 



Af— >oo,A<7 



6472 



both of which can be seen to flt precisely over the whole 
iV-resonant region. 

Having found good approximations for both regimes, 
the last missing piece is the crossover. For even N, the 
smooth shape of the crossover in Fig. 6 suggests a simple 
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FIG. 6: (Color online) A dependence of the contact reflection 
R for various even (upper panel) and odd (lower panel) con- 
tact lengths A*'. Solid: the exact value as given in Eq. (2). 
Dashed: the limit iZjv^oo given in Eq. (4), along with its 
approximation A^/647^, valid for A ^ 1. Dotted: the 
approximation exp (— A'^A/7), valid for even A'' in the A''- 
resonant regime. The inset shows the identical data in the 
semilogarithmic scale, further illustrating the precision of the 
exp {—NA/j) approximation in the TV-resonant regime. 



function of the form R — \/A" 
find that for the case n — 1/2, 
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gives an extremely good match over the full range of A. 
Moreover, a very similar function is found to match the 
crossover for odd values of N, 
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Both approximations show slight deviations from the ex- 
act value for small N but match with high precision for 
larger N . Obviously, the two reflection probabilities be- 
have like squares of quantum mechanical amplitudes in- 



terfering either constructively or destructively with each 
other. 

The parameter values where -Rrosonant ^^nd Rn^oc co- 
incide are of special interest. In the case that N ^ I, 
where this coincidence happens for A ^ 7, the condition 
for this is simply 

exp(-iVA/7) = AV6472, 

leading to an expression for the A-dependent effective 
contact length 



«.,(A) = |..(5 



(6) 



over which a longer contact does not further modify 
transport. This can be interpreted as the length that 
contributes to the electron transmission for a very long 
contact. 

The inverse of Eq. (6) can expressed using the 
Lambert-VF function, ■^^ 



A 



opt(iV) = 2jW{AN)/N 



which can be approximated in the range of interest as 

Aopt(iV) « 27lniV/iV. (7) 

A. Generalization to arbitrary injection energies 

Having found the transport relations at the fixed en- 
ergy E = e, we can now continue with generalizing the 
results for E ^ e. Assuming the general functional form 
of the reflection 



-Rrosonant = CXp (-2iVA/ai) 



we numerically find the following: 

ai - ^4^ - ^2 (1 + C) , 

which holds for arbitrary fixed N with an approximate 
error estimate |CI ^ 1/^ capturing the resonant oscilla- 
tions. Considering the characteristic form of the density 
of states of the linear chain. 



p{E) ^ (n^i^^-E^y 



we can rewrite the last expression as 

ai w I/ttp, 

reflecting the similarity to a weak point contact where the 
tunneling transmission is proportional to the density of 
states on either side. The last relation could be confirmed 
numerically to hold very generally, as will be discussed 
below in the discussion of realistic contacts for nanotubes 
and ribbons. 
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FIG. 7; (Color online) The contact reflection for even A'' ob- 
tained from Eq. (5). Quite visible are the two regimes sep- 
arated by the minimal line Aopt(A''). The "exact" value for 
Aopt is the true minimum for fixed TV. The "approximate" 
value comes from Eq. (7). 
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FIG. 8: (Color online) Generalized model including nondiag- 
onal terms: the individual wideband leads for each atom in 
the extended contact region are replaced by a metal with an 
internal structure, here modeled as a 2D-square lattice. The 
new parameters are 7|| and y±, describing the internal hop- 
ping in the lattice parallel and perpendicular to the contact 
surface, as well as 7c, describing the hopping at the contact. 
For simplicity, we consider only the isotropic case 7|| — 7±. 
This leaves us with the single effective parameter A — 7c /7±. 



IV. NONDIAGONAL CONTACTS 



For the regime of A^ -^- cxi, we can similarly, i.e. 
merically, find the expression 



a2 



7 



(47^ - E^ 



which fits the exact formula [Eq. (2)] with arbitrary pre- 
cision for fixed A ^ 7 small enough and for E not too 
near the band edges. Unlike the formula for ai, how- 
ever, expressing a2 as a function of the density of states 
alone does not help to generalize the relation to other 
structures. 

The crossover region, generally governed by interfer- 
ence effects, can be approximated by averaging over 
quantum mechanical phases, resulting in 



R, 



■crossover 



Rn- 



Rr, 



which gives a good approximation for the full parame- 
ter space with A^ ^ 1, A ^ 7 and E away from band 
edges. Apart from the resonant oscillations, this now al- 
lows the full description of the reflection, and we find a 
precise numerical confirmation of the previously obtained 
expression of the effective contact length. 
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A complete overview of the contact reflection and both 
parameters N and A is shown in Fig. 7. 



To generalize our results beyond the diagonal contact 
approximation, we model the contacting metal not as 
a single-parameter wideband lead but as a material with 
an internal structure, leading to off-diagonal terms in the 
contact matrix (see Fig. 8). Figure 9 illustrates that the 
off-diagonal terms in the self-energy do not bring any 
qualitative changes to the behavior described before. An 
exact quantitative mapping would depend strongly on 
the details of the model. 



V. REALISTIC CONTACTS FOR CARBON 
NANOSTRUCTURES 

For the case of CNTs and GNRs, the method of the 
Chebyshev polynomials cannot be used to obtain an ana- 
lytical solution due to the the nonTrpcommutativity of the 
partial Hamiltonians of the periodic structure. In numer- 
ical studies, however, we find that the behavior is identi- 
cal to that of the linear chain, except for a quantitative 
adjustment of the parameters ai and 02 (see Fig. 10). 
The A^-resonant regime can be described precisely by a 
simple generalization of the law found for the linear chain 

^resonant = exp (-2A^A/ai) , 

where Nch is the number of channels and p the total den- 
sity of states per unit cell. Generally, both values are 
dependent on the energy and the chirality of the tube 
or width and edge geometry of the ribbon. For metallic 
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FIG. 9: (Color online) Nondiagonal contacts: Reflection R 
of the generalized model displayed in Fig. 8. Relating the 
parameter A = 7c /7± to the parameter A of the wideband 
leads, the results are qualitatively similar to those of the orig- 
inal model (Figs. 5 and 6). One prominent difference is the 
enlarged reflection -R(A) for A/7 between 1 and 10: While 
the diagonal self-energy was uniform for every atom along 
the contact, the nondiagonal self-energy now is sensitive to 
the edge of the contact. For large values of A, where only the 
atoms near the edge contribute to the transport, this causes 
the visible deviation from the R — A^/qI law. The same 
reason is behind the visible irregularities in the resonant os- 
cillations of i?(iV). 



CNTs near the Fermi energy, however, one finds the gen- 
eral values of Nch — 2, p ~ 2Nch/'ijdcc, and, therefore, 
ai = 37dcc/27rCc- (7 = 2.66 eV and dec = 1-42 A). 
Introducing the physical length of the contact region 
L = ^uc-^ with the length of the unit cell l^^, the previous 
formula can be rewritten as 
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FIG. 10: (Color online) Extended contacts for a (6,6) CNT 
(left) and the corresponding graphene nanoribbon. Top panel: 
density of states with characteristic van Hove singularities. 
The zigzag-edge state in the nanoribbon causes a peak at 
E — Ep. Center panel: The value ai = — 2A'"A/ln7?, here 
computed for A'^ = 40 and A — 10~^ eV, lies already very near 
to the limiting case ai — Nch/Trp. Bottom panel: The value 
a2 — A/\/i? is well converged for N —* 00 and A = 10"'^ eV. 
Note the suppression of both ai and Q2 in the ribbon at Ep 
where the presence of the localized edge state suppresses the 
conductance in the contact region. 



For the A^-independent regime, the general law of 
Rn^oo — A/q;2 still holds, but the functional form of 
the parameter a2 at arbitrary energies could not be de- 
termined. Generally, it turns out that a2 is suppressed at 
van Hove singularities in a similar way as ai is. Further- 
more, metallic CNTs have a fairly constant value of a2 
around Ep. At E — Ep, we find a2 — 4.247 for armchair 
CNTs and a2 — 5.667 for metallic zigzag CNTs. 

For GNRs, the situation is slightly more complex due 
to the presence of edge states at zigzag edges. ^'*'^° In 
metallic ribbons with armchair edges, the situation is 
similar to that of metallic CNTs, and we find a value 
of a2 — 8.O7 dX E = Ep. The quantitative difference 
from the value of the corresponding CNTs can be ex- 
plained by the presence of only one conduction channel 
at the Fermi energy. For ribbons with zigzag edges, how- 
ever, our simplified model suggests that the constant a2 
is completely suppressed ai E = E-p due to the peak in 
the density of states, caused by the edge state. 

Physically, this suppressed value of both ai and 02 
suggests that the injection of charges into the edge state 
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FIG. 11: (Color online) Differential conductance through a 
two-terminal setup of a linear chain with symmetric, infinite- 
length extended contacts. The central region consists of 
A^conductor = 200 atoHis. The distance between the Fabry- 
Perot oscillations in the gate voltage Vg depends on the total 
gate capacitance Cg as ^Vg = e/Cg. The extent of the di- 
amonds in direction of the bias voltage Vh depends directly 
on the level spacing SE = i'F/i/2Lconductor = 277r/Ai'conductor 
as 5Vh — eSE. Top; strong coupling A — IO7, leading to 
an extremely short effective contact length producing sharp 
resonances and a distinct diamond pattern. Bottom: moder- 
ate coupling A = I7 leading to an effective contact length of 
about four unit cells. The oscillations have a sinoidal shape 
with strongly reduced amplitude. For yet weaker coupling as 
it is to be expected for Pd or Ti, the amplitude is rapidly 
reduced even further making the oscillations undetectable. In 
this case, a stong defect within the contact region may act as 
a point of scattering and recover sharp resonances. 



via extended leads is less efficient than that into the con- 
duction channels of nanotubes or armchair ribbons. How- 
ever, taking into account the results of detailed ab initio 
calculations reveals a spin splitting of the edge states, 
which strongly affects the bands at the Fermi energy and 
opens special spin transport channels,^*' which are not 
captured by our model. 



scattering region. In the zero-bias differential conduc- 
tance, the spacing of the Fabry-Perot resonances depends 
on the gate capacitance Cg alone as 5Vg — e/Cg. Only 
the diamond shapes in a plot of the finite bias differential 
conductance 
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give access to the spacing of the energy levels 5E = 
Hvp/Lq and can thereby be used to measure the length 
Lq of the resonator. 

As visible in Fig. 11, however, the amplitude of these 
oscillations is strongly reduced as soon as the effective 
contact length exceeds the length of one unit cell. One 
could view this situation as smooth contacts that cause 
the Fabry-Perot oscillations to be broadened and the res- 
onator length Lq to be ill defined. 

In some experiments using extended contacts on 
CNTs, the length of the scattering regions was mea- 
sured to be just as long as the uncovered region of the 
tube,^ which could be explained based on our model by a 
strong contact A and, therefore, a short effective contact 
length. For weaker contacts A, it is to be expected that 
Fabry-Perot oscillations cannot be cleanly observed any 
more. A point defect inside the contacted region might, 
of course, act as a scattering point instead and give rise 
to oscillations that indicate a resonator longer than the 
uncovered region of the CNT. 

In our previous study,^^ we chose to average this os- 
cillating conductance over Ep ± 0.5 eV in order to sep- 
arate the finite-length Fabry-Perot effects from the ef- 
fects caused by the contacts themselves. Physically, 
this is similar to the thermal effects caused by high 
enough temperature. For the chosen conductor length of 
Lq — 100 nm, this approach was very successful in can- 
celing all Fabry-Perot oscillations and in reproducing the 
physics of a single extended contact. The resonance os- 
cillations within the contact were, of course, also strongly 
suppressed by the averaging, leaving only a minimal sig- 
nature that we correctly identified as such. 



VII. NONEPITAXIAL CONTACTS 



VI. TWO-TERMINAL SETUP AND 
FABRY-PEROT PHYSICS 

A realistic setup for conduction measurements in CNTs 
and GNRs generally needs a second contact at the other 
end of the system to close a circuit. Such a setup is well 
known to lead to Fabry-Perot-like oscillations of the con- 
ductance along the energy range. ""^ For very bad contacts. 
Coulomb blockade has been observed, but we intention- 
ally avoid this regime that would demand the inclusion 
of charging effects. 

One important aspect of Fabry-Perot oscillations is 
their experimental use in measuring the length of the 



Unlike the theoretical model contacts presented so far, 
realistic samples produced in experiment are never per- 
fectly epitaxial but contain imperfections due to fabrica- 
tion faults, lattice mismatch or metal faceting. To check 
whether the effects described so far are robust to such 
perturbations, we have investigated various kinds of dis- 
order at the contact. A relatively weak disorder was im- 
plemented as random fluctuations of the contact param- 
eter A on each atom i as 



. fluct 



A{l+CiW), 



with an evenly distributed random variable — 1 ^ Ci ^ 
1 and a parameter W specifying the relative strength of 
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FIG. 12: (Color online) Effects of two different kinds of disor- 
der on extended contacts: relatively weak fluctuation disorder 
with varying contact strength A; on each atom i and stronger 
dilution disorder with only a randomly selected fraction of the 
atoms in the contact region attached to a lead. In each case, 
the parameter A refers to the average contact strength. For 
details, see the text. 



the fluctuations. As can be seen in Fig. 12, even for the 
strongest possible value W = 1, the effect of the disorder 
is moderate and purely quantitative. 

An even stronger disorder was realized by using a 
model of diluted contacts, where only a randomly se- 
lected fraction of the atoms in the contact region is con- 
tacted. 



L diluted 



A/P with probability P 
with probability 1 — P. 



This kind of disorder modifies to the observed behavior 
to a much larger degree, but even in the extreme case of 
a 1% dilution (i.e. P = 0.01), the general trend of the 
original model is well preserved (see Fig. 12). 



VIII. MATERIAL RELATED CALCULATIONS 

To link the model results obtained so far to the physical 
properties of real contact materials, we performed den- 
sity functional theory (DFT) calculations of Ti and Pd 
monolayers interacting with a graphene layer as described 
before.^*' We described the valence electrons by TrouUier- 
Martins pseudopotentials and used the Perdew-Zunger 
form of the exchange-correlation functional in the local 
density approximation to DFT, as implemented in the 
SIESTA code.''''^ With a double-zeta basis and a 100 Ry 
energy cutoff in the plane-wave expansions of the elec- 
tron density and potential, we found the total energy to 
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FIG. 13: (Color online) Analysis of the hybridization between 
a graphene sheet and a metal monolayer. As visible in the 
upper scheme, the hybrid band structure matches well with 
an overlay of the band structures of the two individual sys- 
tems rigidly shifted in energy and hybridization at some band 
crossings. Highlighted in the Pd/C and Ti/C band structures 
are the regions of interest, i.e., those hybridizations that con- 
tribute most to the electron injection. 



be converged to < ImeV/atom. We performed a full 
structure optimization to determine the equilibrium ad- 
sorption geometry, the adsorption energy, and the local 
charge redistribution caused by the metal-graphene in- 
teraction. Since the interatomic distances in bulk Pd 
(2.7 A) and Ti (2.95 A) lie close to the honeycomb spac- 
ing in graphene (2.46 A), we considered only epitaxial 
adsorption. For both Pd and Ti, we found a slight pref- 
erence for the sixfold hollow site on graphite. For Pd, 
we found the equilibrium interlayer distance to be 3.2 A, 
consistent with a relatively weak, mostly covalent bond 
energy of 0.3 eV per Pd atom. The interaction between 
an epitaxial Ti monolayer and graphene was only insignif- 
icantly stronger with 0.4 eV per Ti atom at an interlayer 
distance of 3.0 A. 

To study the electronic coupling between the two sys- 
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terns, we first inspected the band structure (see Fig. 13). 
Especially for Pd as a contacting metal, the extraction 
of parameters for our model is greatly simplified by the 
fact that the band structure of hybrid lies close to the 
superposition of the metal and carbon band structures. 
One can see a rigid shift of the carbon bands by Ec — 
0.374 eV while the palladium bands are shifted slightly in 
the opposite direction with AEp^ = —0.020 eV. On top 
of this rigid shift, one can observe slight hybridization 
effects in the band structure. For injecting conduction 
electrons into a graphene sheet or the wall of a carbon 
nanotube, the most important area of the Brillouin zone 
is the K point, the Fermi point of graphene. In Fig. 13, 
a small avoided crossing is visible near this region in the 
Pd/C band structure. To extract an estimate of tight- 
binding parameters from these data, we modeled a hon- 
eycomb lattice and a matching hexagonal lattice repre- 
senting both sheets. As it turned out, a single orbital per 
atom is sufficient to obtain bands that can be fitted to the 
hybridizing bands near the Fermi level with a single pa- 
rameter each. Now, an additional coupling between the 
two sheets was introduced, linking each Pd atom with its 
six neighboring C atoms. This hopping parameter could 
then be tuned to reproduce a hybridization between the 
two subsystems, which is close to that in the hybrid band- 
structure obtained from DFT, resulting in a coupling of 
ipd/c « 0.15 eV. 

For the case of Ti, the distortions in the band structure 
caused by the hybridization of the two layers are consid- 
erably stronger than for Pd. Still, the change in the 
carbon related levels can be modeled by a rigid shift of 
AEq — —1.15 eV. To determine the hopping parameter, 
the same procedure as for Pd could not be directly ap- 
plied, because the relevant band of the Ti-monolayer can- 
not be reproduced with a single-orbital hexagonal lattice. 
Instead, a rough estimate was obtained by visually com- 
paring the bandstructures themselves where the avoided 
crossing near the -RT-point is at least twice as large as for 
Pd, giving an estimated value of txi/c ^ 0-3 eV. 

To turn these parameters into values of A that can be 
directly placed into our model calculations, we need the 
surface density of states, which is comparable for both 
materials at A/mo ~ 1 eV~^. Finally, the connectivity 
at the interface is also important: Each C atom con- 
tacted to three different metal atoms simply triples the 
value of A. The internal connections inside the metal 
are already taken into account with the surface density 
of states and do not have to be considered any further. 
With the relation A — t^M , this again gives rough esti- 
mates of Apd ~ 0.06 eV and Ati > 0.3 eV. 

Unfortunately, this approach of computing a graphene 
layer and a layer of the contacting material within a com- 
mon unit cell cannot necessarily be transferred to other 
materials of interest in any straightforward way. There 
exist, however, ah initio calculations of various metals 
in contact with graphene or CNTs that show a clear 
trend: ^^'■^^'^^ The highly conducting metals Au, Ag and 
Cu generally have a very weak binding energy, insuffi- 



cient for wetting the carbon surface; so, a clean contact 
is hard to achieve. Pd, Pt and Ti all have sufficient bind- 
ing energies for wetting the surface. Pt and Ti both have 
higher binding energies than Pd. For pointlike end con- 
tacts, such strong bonds give good transparency.^'' For 
extended contacts, however, the weak bonds of Pd arc to 
be preferred. 

Previous calculations^^ attributed the difference be- 
tween Ti and Pd metal contacts to the formation of 
different-sized metal clusters at the interface. Our re- 
sults, presented above, offer a more fundamental ex- 
planation: It is exactly the weak bonding between Pd 
and graphene or CNTs — ^just large enough to wet the 
surface — that makes Pd such an excellent contact mate- 
rial. 



IX. CONCLUSIONS 

To conclude, we introduced a model for electrical con- 
tacts for carbon nanotubes and graphene nanoribbons 
that captures the fact that the contacts in an experi- 
mental setup typically extend over a length of several 
tens of nanometers, covering the carbon structure with 
some contact metal. We have demonstrated the coun- 
terintuitive result that, given a metal coating of several 
nanometer length, the contact transparency is actually 
improved by using a metal that couples more weakly to 
the carbon surface. Using ab initio results of Ti and Pd 
as contact metals, we have demonstrated that Pd actu- 
ally forms a weaker bond, giving an explanation of the 
experimental finding that Pd forms good contacts. This 
finding suggests a possible route for future attempts in 
optimizing the charge injection in carbon nanotubes and 
graphene, that is, to find contacting materials that cou- 
ple to the carbon surface as softly as possible to exploit 
the available contact length. 

Starting with a detailed analysis of an analytically 
solvable minimal model, we have demonstrated by nu- 
merical calculations that the qualitative results are ro- 
bust to various modifications of the system. Replacing 
the atomic wire by the actual atomic structure of a car- 
bon nanotube or a graphene ribbon can be accomodated 
by adjusting just two parameters at any given energy. In- 
cluding a second contact to model a realistic two-terminal 
conductance measurement leads to Fabry-Perot oscilla- 
tions that can be averaged out to give the original result. 
A weak disorder in the contact interface has very little 
quantitative effect, and even a strong disorder (similar 
to metal grains forming contact only in certain points) 
leaves the qualitative behavior unchanged. 

In view of a possible experimental confirmation of 
our results, the biggest challenge may lie in the fabrica- 
tion of well-controlled finite-length contacts down below 
the magnitude of the effective contact length of a few 
nanometers. With current technology, this precision is 
yet out of reach, but with future developments, it should 
well be possible to tune the geometry of contact with 
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sufficient precision. An alternative approach would be to 
tune the coupling contact interface by some means. Di- 
rect tuning of the bonds, as it can be done for molecular 
junctions, '^'^ seems unfeasible for this kind of geometry. 
Instead, the insertion of an insulating atomic layer below 
the contacting metal is already being used to improve the 
contacts for graphene monolayers and might work as well 
for nanotubes. 
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APPENDIX A: ANALYTICS FOR THE 
ONE-DIMENSIONAL MODEL 

1. Transmission calculations 

The Hamiltonian of a two-probe system for transport 
calculations is given by 



with 



n 



i/L i^Lc 

H^c Hji 



(Al) 



where He describes the finite-size conductor region and 
Hl/r describes the leads, which are connected to inde- 
pendent reservoirs and have no direct contact with each 
other. From 7i = TiMt follows that H^^ — H\^ and 

^Rc = ^cR- 

To simplify the notation, we first define the complex- 
energy Green function: 

g{£) = {£-ny^ 

and derive from it the expressions for the retarded and 
advanced Green functions {£ = E ± irj) 

g'{E) ^ lim g{E + iT]), 
g^E) = lim g{E-irj). 

The transmission through this system is given by^'^'^"'^ 

T = TY{T^giTngn, (A2) 
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The Hamiltonian of the model at hand, depicted in Fig 1, 
can be split up according to Eq. (Al): The conductor 
consists of just one atom, so its Hamiltonian is a 1 x 1 
matrix H^ = ( e ) . The left lead is a semi-infinite chain, 
contacted only at the last atom. 



Ht = 



H, 



cL 
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The right lead consists of a chain of A^ atoms, each at- 
tached to a wideband lead. This can be captured by 
defining an effective Hamiltonian of the form 



rrcff _ 
^R - 
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(A3) 
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together with a contact point in the first atom only, 
^cR = (-7 ... 0)^^^. 

Note that H^ is the effective Hamiltonian containing 
the retarded self-energy, so g\^ = {E ^ iO+ — H^^ and 

t;g = (s-io+-(iJf)t)-\^ 

In the following, we will simplify the notation by set- 
ting 7=1 and £ = 0. Both constants can be rein- 
troduced in the final result [Eq. (A8)] by substituting 

E^{E-e) h. 

2. Inverse based on Chebyshev polynomials 

As a starting point for a full analytical solution, we 
look at a finite linear chain of length N, which has the 
Hamiltonian 



H^ = 



/ -1 
-1 

-1 
V ... 



\ 



-1 

-1 
-1 / 



NxN 



The quantity of interest of this system is the 1 , 1 matrix 
element of the Green function g^ {£) = {£ — Hn)~ ■ 
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The solution is based on the Chebyshev polynomials of 
the second kind,^^ which can be defined via the determi- 
nant identity, 
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or, equivalently, by the recursive definition, 

Uo{x) = 1, 
Ui{x) = 2x, 

Un+l (x) = 2xUn (x) ~ C/„_i (x) . 



(A4) 



We can now use the well-known identity for the matrix 
inverse, 
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3. Surface of semi-infinite linear chain 

The surface Green function of a semi-infinite linear 
chain can be defined as 



gs{£) = lim [g"" {£)],, 

= lim fN{£/2) 

=: /oo(f/2). 

To find an expression for fo^ (x) , we can use the recursive 
definition of the Chebyshev polynomials [Eq. (A4)] and 



For N ■ 



In (x) = {2x - JN-i {x)) ^ . 
cx), this becomes 

foo{x) = (2x-/oo (x))"\ 

which has two solutions foo {x) = a; I 1 ± -^/l — 1/a;^ 

On the real axis, it follows from Eqn. (A4) by induction 
that l/jv (a;) I < 1 when |a;| > 1, so that we can select the 
correct solution 



/oo(a;) 



1 - Vl - 1/a 



(A6) 



which can be continued analytically to x £ C\(— 1,1) 
by reading the square root of a complex number as the 
principal square root, uniquely defined everywhere except 
on the negative real axis by the condition Re {\/x) ^ 
0, Va; G C. 

The retarded surface Green function follows as 



g^ {E) = hm gs {E + iry) 

r;— ►0+ 



I E/2- ^^2/4-1 for \E\ ^ 2 
E/2 - i^l -^2/4 for \E\ < 2 



(A7) 



4. Transmission of the model system 

With these results, we can now obtain the quantum 
mechanical transmission of our model system. The left 
lead is a semi-infinite chain giving a self self-energy of 



i?cL(f 



Hv) ^Ht 



For |£:| ^ 2, Gl [E) is real [see Eq. (A7)], so, Tl and with 
it by Eq. (A2) the whole transmission T, are strictly zero. 
In the following, we therefore assume |i?| < 2 and select 
the second case in Eq. (A7): 



Sf. = 



E 



1-1^4/^^-1 



To find the self-energy of the right lead, we use the defi- 
nition of H"^ from Eq. (A3) and find 

-1 



^K 



= H,n{E~H''+iA/2) 
= /jv (S/2 + iA/4) . 
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Now, we can put together all parts to calculate the trans- 
mission, 



5! = 

Tl = 

Tr = 

T = 



2(e + iE^^A/E^ - 1 - 2/jv {E/2 + iA/4) 



2(E- lE^AjE'^ - 1 - 2/jv (-B/2 - iA/4) 

E^jA/E^ - 1, 

i [/at {EI2 + iA/4) - Jn {EI2 - iA/4)] 

2Im[/^(£;/2-iA/4)] 

8V4 - E"^ Im (/at (EI2 - iA/4)) 
I E - iV4 - ^2 _ 2fN (-B/2 - iA/4) I ^ 
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